The number of capillary bridges in a wet granular medium 
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We observed the appearance of capillary bridges in a granular medium consisting of glass beads 
after adding small amounts of liquid. We found the initial bridge formation depending on the bead 
roughness. Furthermore we obtained a statistics for the average number of bridges for randomly 
packed beads in dependence of the liquid content and were able to find an explanation therefore 
based on recent models and former experimental data. 

PACS numbers: 05.65.+b,45.70.-n,45.70.Mg 



I 

o 
o 



> 

O 
m 
o 



X3 
O 

o 



X 
J3 



If one adds a wetting liquid to a granulate, there oc- 
cur capillary bridges which exert forces on the particles. 
These forces are the main reason for the strongly changed 
mechanical properties of wet granular matter compared 
to the dry case 0, S H S H [1 ■ Despite the notable 
advance in the comprehension of the dynamical prop- 
erties of dry granular materials 0, ^] , the physical 
mechanisms, which are the cause for the character of wet 
systems, remain mostly unknown. 

Since the formation of the liquid bridges depends ap- 
parently strongly on the microscopic geometry of a gran- 
ulate, an adequate description of the packing is essen- 
tial. The random close packing of beads, which is the 
simplest random packing and which occured in our ex- 
periments presented in this paper, is an interesting and 
challenging topic in itself and has been much studied in 
the past [nlllllll. It was the goal of our experiments 
to show how the formation of capillary bridges, which 
are essential for the mechanics of a wet granular mat- 
ter 

BIS II 111 El, depends on the geometrical particle 
distribution, which was subject of earlier studies .12,, 13] . 

In order to observe and count the single bridges we 
used an index matching technique ^S|- At this we put 
glass beads with diameters of 375^m or 555/Ltm, which 
had a refraction index of 1.51, into a mixture of toluene 
and diiodomethane with the same refractive index, so 
that the granular matter became quite transparent and 
we were able to zoom with the microscope through the 
glass bead layers and observe the bridges. In order to 
achieve this index matching the volume share of the 
toluene was 88.1% and of the diiodomethane 11.9%, re- 
spectively. This liquid mixture took over the role of " air" 
in a "normal" wet granular matter. As the wetting liq- 
uid, i.e. as liquid which forms the capillary bridges, we 
added small amounts of water with fluorescein. 

The granular matter was put into a cuvette (0.95 x 
0.95 X 4cTO, filling height: 2.5 — 2>cm) and was shaken in 
small horizontal circles with an amplitude of 5mm and 
with a frequency of roughly 30-ffz for some minutes, until 
the "wetting liquid", clearly visible due to fluorescein, 
appeared homogeneously distributed. The diameter of 
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the glass beads used was 375/im and 555/i77i. In order 
to prevent crystallization, the small beads were chosen 
to be slightly polydisperse: the spread in bead size was 
10% for the glass beads with a diameter of 375/im and 1% 
[2lj | for the glass beads with a diameter of 555/im. After 
shaking, the volume fraction of the glass was 0.57, i.e. we 
found the granulate packed less densely as it would be the 
case for the loose random packing of 0.60 for the dry case 
[Tll | . We could also reach a random dense packing of 0.62, 
which is the other extreme of the random packing of glass 
beads and is similar to a volume fraction of the glass of 
0.64 =: prdp E3> by tapping the sample several times. 

After adding very small amounts of liquid, there was 
no formation of capillary bridges. The reason for this be- 
havior is that the liquid was first trapped on the surface, 
as it can be seen in Fig. la. Only at a water content of 
0.07%, bridges on almost every bead contact were formed 
(see below). Fig. lb illustrates the case for fully formed 
bridges. 

It is well known that bridges need a certain amount of 
liquid content W to form properly, since the liquid is at 
first bound on the bead surfaces due to the roughness 0, 
0, 01 . The connection between the roughness amplitude 
5 and the saturation bridge volumell^ is given by Wh ~ 
2R5^, where R is the bead radius, and Wf, the volume of 
a bridge at the water content of saturation Wb Our 
measurements |l8j suggest that the relationship between 
the water content W and the volume of a single bridge 
w is generally given hy w = ar^W, where a w 0.25. Q 

From Wb ~ 0.07% we get i5 ^ 500nm, which is similar 
to the peak-to-peak roughness we found from the inspec- 
tion of the beads by atomic force microscopy. 

We obtained the average number of bridges per bead 
by zooming with the microscope through the sample and 
counting the bridges. We carried this out for different 
liquid contents for the random loose packing which was 
in our case 0.57 as well as for the random dense packing 
of 0.62. In order to reach better statistics we counted the 
bridges of 40 beads for each data point for the dense pack- 
ing and 20 beads for the loose packing and calculated the 
average bridge number N. We took the standard devia- 
tion as error for N. Since we counted only bridges from 
beads of the second to the sixth layer, an estimation of 
the influence of the sample boundary is necessary. There- 
fore we measured N for different layers L. We found that 
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Glass bead diameter 0.375 mm 
□ random loose packing 
■ random dense packing 

Glass bead diameter 0.555 mm 
o random loose packing 
• random dense packing 
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FIG. 2: Average number of capillary bridges A'' in dependence 
of the water content W 



FIG. 1: Fluorescence microscope pictures for different water 

contents W; glass bead diameter: 375/im 

la) W — 0.055%, no real bridges have formed yet 

lb) W — 0.2%, capillary bridges between the beads are clearly 

visible 

N{L) = const for 3 < i < 6, whereas for L < 2 sniall 
variations were observed. We believe therefore that we 
can exclude any effect of the sample boundary on TV. 

The dependence of N on the water content is shown in 
Fig. 2. N is equal to zero for very small W. At a water 
content slightly smaller than Wb it goes up very fast, then 
the curve becomes flatter. It can be seen further, that 
N{W) for the random dense packing is about 10% larger 
than N(W) for the random loose packing. 

As it was mentioned earlier, the number of capillary 
bridges is influenced strongly by the microscopic geomet- 
rical properties of the granular medium. A model which 
deals with the number of contacts of randomly packed 
glass beads is the caging model Tsj. According to this 
model one needs 4.79 =: Ndc bead contacts on average 
in order to pin a bead with other beads. We assumed 
that this model provides at least roughly an appropriate 
explanation for our experiments and have defined Wb ac- 
cordingly. Consequently Wb is the water content at which 
a capillary bridge exists then and only then, if two beads 
touch, i.e. if their surface distance d is zero. The smaller 
non-zero values a.t W < Wb can be explained as follows. 
The water film on the beads varies in thickness at differ- 



ent areas of the bead surface. Therefore it is possible at 
certain water contents that bridge formation occurs only 
at a fraction of the bead contacts. On the other hand, for 
W > Wb there can also occur slightly " streched" bridges 
between beads which do not touch, but whitch are close 
enough to each other (see below). 

Our assumption that the caging model is appropriate 
for our case is confirmed further if we compare the mea- 
sured distribution of the capillary bridges at water con- 
tents W « Wb for the random loose packing to the distri- 
bution obtained by a simulation 13] based on the caging 
model. As it can be seen in Fig. 3, the measured curve is 
very similar to the simulated one. Particularly the mea- 
sured curve does not exceed the value of 8, i.e. it fulfills 
a necessary condition for the validity of the caging model 
[l3l |. At higher water contents the histogram curve is 
shifted to the right, but we did never see a glass bead 
with more than 9 capillary bridges. |23j. 

We will compare now our measurements concerning the 
dependence of the number of capillary bridges on W to 
the measured distance distribution of beads in a random 
close packing |0. For W > Wb there exists at least at 
every bead contact (i.e. at every point with c? = 0, where 
d is the surface distance of the beads) a capillary bridge: 
N > Nc- Nc is here the number of real contacts per 
bead. According to the caging model we have N^c = 4.79. 
[l3l | Bridges can also exist between beads which do not 
touch, but which are only sufficient close to each other. 
Due to the hysteretic nature of bridge formation ^ , it 
is reasonable to suppose that such bridges can form only 
between beads which collided in the past and then moved 
away from each other. For sufficient strong shaking this 
should be statistically the case for a certain part A* of 
the beads, i.e. 

N{W) = Nc + A*Nb{W) 
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FIG. 3: Distribution of the number of bridges: 

measured (W — 0.15%, random loose packing; open symbols) 

simulated (closed symbols) 



Here Nb is the number of nearest neighbors, which do not 
touch the bead, but whose distance d is smaller than the 
bridge rupture distance dp. The total number of neigh- 
bors in between a distance dp is 



N{W) - N, 
A* 



Nt{W) = N, + Nk{W)=N, 

Herewith we get 

N(W) ^A*Nt{W) + {l- A*)Nc 



(1) 



We need now to correlate the water content W to rup- 
ture distance dp. According to this correlation is 
approximately given by 

1 1 W3 

dn — H 

^ 10 i? 

where R is the bead radius and w the b ridg e volume. 
In accordance to our own measurements |24| . we have 
w — aR^W with a « i. Therefore 



dp = R[{aWy^ + —{aW) 



Now it is possible to calculate N from the measurements 
of A^t 12], plot it over dp and compare it with the directly 
measured number of capillary bridges. Before doing this 
we defined W* := W — Wb, i.e. the water content above 
Wb which is located in the bridges. The reason therefore 
is that the glass bead surface absorbs a certain amount of 
liquid Wb before bridge formation starts. Accordingly we 
defined d* := dp(W*) as corresponding rupture distance. 

In order to get a quantitative expression the values of 
A* one has to consider the mechanism of bridge forma- 
tion. A bridge between two non-touching beads can only 



FIG. 4: Comparison of the measured number of bridges A'^ 
with random loose and random dense packings to the A'' which 
were calculated from the literature values \\^ according Eq. 
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Inset: Number of bead contacts for different packing densities. 



exist if their distance is smaller than dp if this beads 
touched each other in the past. Therefore one would 
expect A* = J for a one-dimensional system with statis- 
tically moving particles. In a real granulate the situation 
is more complicated. It is possible that particles move 
against each other by a distance smaller than dp without 
ever touching, which leads to a decrease of the value of 
A*. 

In order to obtain an estimation for A* we consid- 
ered a single bead with radius R in the point of ori- 
gin. Then we assumed other beads with radius R mov- 
ing by at an impact parameter p* which was chosen 
at random. It is evident that bridges can only form if 
< p* < 2R, i.e. a bead collision is a necessary con- 
dition for bridge formation. If the moving bead moves 
along the x-axis, the distance between the surfaces is 
given by c? = \/ x*^ + p*^ — 2R. At a certain point of 
time a bridge can only exist if d < d*. A further condi- 
tion for the existence of a bead is that the collision took 
back in the past which means that only half of the beads 
with < p* < 2_R and d < d* will have a bridge. The 
arc length of a bead which has a bridge is given by 



Blip* 




with xi = v/(2i?)2 and X2 = ^ (d* + 2R)^ -p*^. 

The arc leangth of a (bridgeless) bead with 2R < p* < 
2R + dp and d* < d* is analogous 



Blip*) = 2 
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With h 



2R+dl 
2R 



:= 2np* J^'' Bi{p*)dp* and h 
B2{p*)dp* one finally gets 

h 



A*id;) 



2/1 + 12 



A* (dp) can be approximated very well as follows: 

1 dl 
~ 1.447— + 2 

A{d;/R) R 

Fig. 4 shows the plots of our measured N over d* for the 
beads with a diameter of 555/im and for the random loose 
packing and the random dense packing. The drawn curve 
represents the values of N for the dense packing which 
were calculated from the literature values of Nt |l2l | using 
Eq. n It is clearly visible that our calculated values lie 
near to those of the literature for Nc = 3.3 for the loose 
and Nc = 4.3 for the dense packing. Nc differs from 



the Nc of the caging model significantly for the loose 
packing, while the difference in the dense packing is quite 
small. Hence the assumption that Nc depends directly 
on the packing density pp is self-evident. The measured 
Nc{pp) as well as the theory values are in agreement with 
the linear curve N^c — Nc{pp) = 22{prdp — Pp) (see Fig. 
4 (inset)). This plot displays the Nc of the loose as well 
as the dense packing for the used beads with a diame- 
ter of 555 and 375/^to. The point of origin corresponds 
to the Ndc of the caging model for the random dense 
packing prdp- Concluding the value of the number of 
contacts according to the caging modell N^c is given 
by the linear extrapolation of the Nc{pp) obtained from 
our measurements: Ndc — Nc{prdp)- 
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